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NONLINEAR  EVOLUTION  OF  THE  KELVIN-HELMHOLTZ  INSTABILITY 
IN  THE  HIGH  LATITUDE  IONOSPHERE 


I .  INTRODUCTION 

The  Kelvin-Helmholtz  instability  has  played  an  important  role  in 
understanding  a  variety  of  space  and  astrophysical  phenomena  such  as  the 
dynamics  of  comet  tails  [Ershkovich,  1976;  Brandt  and  Mendis,  1979],  the 
interaction  between  adjacent  streams  of  different  velocities  in  the  solar 
wind  [Parker,  1963;  Sturrock  and  Hartle,  1966;  Jokipii  and  Davis,  1969], 
low  frequency  fluctuations  in  the  polar  cusps  [D'Angelo,  1973;  D'Angelo  et 
al.,  1974],  auroral  arc  dynamics  [Hallinan  and  Davis,  1970],  and  the 
interface  between  the  solar  wind  and  magnetospher ic  boundary  [Dungey, 
1955;  Parker  1958,  Aubry  et  al.,  1971;  Lepping  and  Burlaga,  1979; 
Southwood,  1968].  Hones  et  al.  [1981]  have  detected,  using  ISEE  satellite 
data,  vortical  structures  in  the  low  latitude  boundary  layer.  They  show 
that  the  helicity  and  approximate  spatial  extent  are  consistent  with  the 
Kelvin-Helmholtz  instability.  Williams  [1980]  and  Sckopke  et  al.  [ 1 98 1 ] 
have  observed  oscillations  of  the  dayside  magnetopause  with  time  scales  on 
the  order  of  several  minutes  which  are  consistent  with  Kelvin-Helmholtz 
driven  waves. 

Recently,  much  experimental  [Basu  et  al.,  1987;  Weber  and  Buchan, 
1081;  Bythrow  et  al.,  1984;  Cerisier  ct  al.,  1085;  Rodriguez  and 
Szuszczewicz,  1984;  Curtis  et  al.,  1982;  Bake;  et  al.,  1986;  Vickrey  et 
al . ,  1980]  and  theoretical  [for  recent  reviews,  see  Keskinen  and  Ossakow, 
198.5;  and  Kintner  and  Scyler,  1085  and  refetences  therein]  attention  has 
been  given  to  the  origin  of  high  latitude  ionosphet  i r  and  magnet osphci i c 
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plasma  turbulence.  The  Kelvin-Helmholtz  or  velocity-shear  driven 
instability  can  lead  to  both  electric  field  and  density  fluctuations  in 
the  high  latitude  near  earth  space  plasma  [see,  e.g.,  Kintner  and  Seyler, 
1985].  Studies  of  velocity-sheared  flows  in  space  plasmas  can  be  divided 
into  tvo  groups  depending  upon  vhether  plasma  flov  velocities  are  either 
parallel  [Lerche,  1966;  Southvood,  1968;  D'Angelo,  1973;  Ershkovich,  1979; 
Gary  and  Schwartz,  1980;  Buba,  1981;  Mishin,  1981;  Lee  et  al.,  1981; 
Walker,  1981;  Keskinen  and  Huba,  1983]  or  perpendicular  [Hallinan  and 
Davis,  1970;  Miura  and  Sato,  1978;  Miura  and  Pritchett,  1982;  Pritchett 
and  Coroniti,  1984;  Thompson,  1983]  to  the  ambient  magnetic  field.  Both 
cases  have  been  studied  in  the  MHD  [Mikhailovskii ,  1974;  Sen,  1964; 
Southvood,  1968]  and  electrostatic  [D'Angelo,  1965;  Smith  and  von  Goeler, 
1968]  limits.  Furthermore,  the  velocity,  in  both  cases,  is  usually  taken 
to  vary  spatially  transverse  to  the  magnetic  field  in  the  electrostatic 
limit.  In  this  study  we  restrict  ourselves  to  sheared  flovs  perpendicular 
to  the  geomagnetic  field.  Hallinan  and  Davis  [1970]  and  Webster  and 
Hallinan  [1973]  have  attributed  the  small  scale  vortex  configurations 
often  seen  near  auroral  arcs  [Hallinan  and  Davis,  1970;  Oguti,  1974]  to  be 
driven  by  a  transverse  Kelvin-Helmholtz  or  velocity  shear  driven 
instability.  Kintner  [1976]  and  Kelley  and  Carlson  [1977]  have  shown  that 
electrostatic  wave  activity  is  well  correlated  with  strong  transverse 
velocity  shears  associated  with  auroral  arcs.  Recently,  Basu  et  al. 
[1987],  using  Dynamics  Explorer  satellite  data,  have  measured  the  spectrum 
of  the  electric  field  and  density  fluctuations  associated  with  velocity 
sheared  regions  near  auroral  arcs. 

Theoretical  analyses  of  the  transverse  Kelvin-Helmholtz  instability 
in  space  plasmas  have  exploited  the  analogy  with  crossed-f ield  electron 
beams  in  laboratory  plasmas  [Buneman  et  al. .  1966;  Levy  and  Hockney,  1 968 ] 
and  velocity  sheared  flows  in  neutral  fluids  [Chandrasekhar,  1961; 
Michalke,  1964;  Blumen  et  al.,  1975;  Drazin  and  Davey,  1977  ) .  Miura  and 
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Sato  11978]  have  studied  the  transverse  Kelvin-Helmholtz  instability, 
using  numerical  simulations,  and  its  role  in  auroral  arc  deformation. 

Thompson  1 1983 ]  and  Satyanarayana  et  al.  f 1 987a ]  have  investigated  the 
effects  of  parallel  electron  dynamics  on  the  transverse  Kelvin-Helmholtz 

instability.  Sydora  et  al.  [ 1984 J  have  analyzed  the  possible  role  of  the 
transverse  Kelvin-Helmholtz  instability  in  radially  injected  plasma 
clouds.  Pritchett  and  Coroniti  (1984]  have  discussed  the  collisionless 
transverse  electrostatic  Kelvin-Helmholtz  instability  with  inclusion  of 
lowest  order  finite  Larmor  radius  effects.  Recently,  Lyons  and 
Valterscheid  [1985]  discussed  the  generation  of  auroral  omega  bands  in 
terms  of  a  Kelvin-Helmholtz  instability  in  thermospheric  neutral  wind 
patterns. 

However,  the  previous  theoretical  studies  have  treated  the  Kelvin- 
Helmholtz  instability  as  it  occurs  in  a  collisionless  magnetospheric 
context  with  no  collisional  ionospheric  coupling.  Due  to  the  relatively 
long  growth  times  of  the  Kelvin-Helmholtz  instability  compared  to  the 

f 

Alfven  travel  times  between  the  inner  magnetosphere  and  ionosphere,  and/or 
long  static  electric  field  mapping  distances,  ionospheric  collisional 
effects  must  be  considered  in  a  proper  description  of  the  Kelvin-Helmholtz 
instability  in  the  high  latitude  near-earth  space  plasma. 

In  this  report,  we  present  the  first  numerical  simulation  of  the 
nonlinear  evolution  of  the  electrostatic  Kelvin-Helmholtz  instability  with 
ionospheric  Pedersen  conductivity  coupling.  We  find  that  the  Kelvin- 
Helmholtz  instability  develops  in  a  distinctly  different  manner  in  the 
nonlinear  regime  with  Pedersen  coupling  than  without  it.  Specifically,  we 
show  that  Pedersen  coupling  effects,  in  conjunction  with  a  neutral  wind 
and  density  gradient,  (1)  result  in  an  increased  time  scale  for  Kelvin- 
Helmholtz  instability  wave  growth.  < .' )  inhibit  Kelvin  Helmholtz  vortex 
formation,  (3)  lead  to  nonlinear  s  r  •  u-- »  mi  ^  -  •  h  i  h  '.an  no  described  as 
"breaking  waves",  and  (4)  generate.  in  »h*  I'liiimai  i  -v. :  ••••,  small  scale 
turbulence  by  means  of  secondary  instabilities  growing  on  primary  waves. 


Ve  have  also  computed  the  spatial  power  spectra  of  the  electrostatic 

potential  and  density  fluctuations  and  find  differences  with  and  without 

Pedersen  conductivity  effects.  For  the  standard  Kelvin-Helmholtz 

instability  with  a  velocity-shear  gradient  in  the  y-direction,  but  with  a 

density  gradient  and  no  ion  collisions,  we  find  P(kx)  *  kx~3'^  *  and 

P(ky)  «  k  “  *  ±  ‘  for  the  electrostatic  potential  fluctuations,  and 

P(k  )  «  k  ~2-2  *  0,3  and  P(k  )  «  k  _2,3  *  °'3  for  the  density 

xx  y  y 

fluctuations.  Here  P(k  )  and  P(k  )  are  the  transverse,  time  and  space 

x  y 

averaged  spatial  power  spectra.  For  the  Kelvin-Helmholtz  instability  with 
a  density  gradient  and  including  ion  collisional  effects  we  find,  in 
general,  shallower  spectra  with  P(k  )  *  k  ~  5.0)^  ^  j  «  ^  - 

xx  y  y 

for  the  potential  fluctuations,  and  with  P(kx)  *  kx_^‘^ 
2’3^  and  P(ky)  «  ~  2*^  for  the  density  fluctuations. 

The  organization  of  the  paper  is  as  follows.  In  Section  2  ve  present 
the  basic  assumptions  of  the  model,  derive  the  basic  model  equations,  and 
possible  equilibrium  configurations.  In  Section  3  we  present  a  general 
linear  theory  of  the  Kelvin-Helmholtz  instability  with  and  without 
Pedersen  coupling  effects.  In  Section  4  we  solve  the  fundamental 
nonlinear  equations,  in  real  configuration  space,  describing  the  Kelvin- 
Helmholtz  instability  with  Pedersen  coupling.  In  Section  5  we  discuss  the 
evolution,  in  k-space,  of  the  Kelvin-Helmholtz  instability  with  and 
without  Pedersen  effects.  Finally,  in  Section  6  we  summarize  our  findings 
and  apply  our  results  to  velocity  sheared  configurations  in  space  plasmas. 
II.  MODEL  EQUATIONS  AND  EQUILIBRIA 

We  model  the  evolution  of  the  Kelvin-Helmholtz  instability  transverse 
to  the  ambient  geomagnetic  field  and  inrlude  lowest  order  effects  of 
magnetosphere- ionosphere  coupling.  :-inc«  v<=  use  M.u  electrostatic 

approximation,  it  will  be  assumed  that  the  time  scales  of  interest  in  the 


problem  are  long  compared  to  the  transit  time  of  an  Alfven  wave  across  the 
modeled  region  along  the  magnetic  field.  This  assumption  implies  that  the 
plasma  is  quasi-neutral  everywhere  and  that  the  electric  field 
perpendicular  to  the  magnetic  field  maps  perfectly  along  the  field  within 
the  modeled  region.  Further,  the  fact  that  the  electrons  have  a  greater 
mobility  parallel  to  the  magnetic  field  than  the  ions  allows  the 
assumption  that  all  parallel  currents  are  electron  currents  and  that  the 
ions  have  no  velocity  parallel  to  the  field.  As  a  result  of  these 
assumptions,  the  model  may  be  viewed  as  a  set  of  ion  layers  perpendicular 
to  the  magnetic  field  representing  both  ionospheric  and  magnetospheric 
contributions,  as  shown  in  Fig.  1,  which  are  strongly  linked  by  the 
mapping  of  the  perpendicular  electric  field  between  layers  and  by  the 
parallel  electron  currents  which  flow  between  layers  in  order  to  preserve 
quasi-neutrality.  For  simplicity,  we  will  neglect  any  geometric  effects 
due  to  field  line  curvature  and  assume  that  the  magnetic  field  B  =  B  ez  is 
in  the  z-direction  and  constant,  so  that  each  ion  layer  is  an  xy-plane. 
This  assumption  does  not  prohibit  the  magnetic  field  from  taking  on 
different  magnitudes  on  different  layers,  since  it  is  possible  to  show 
that  the  geometry  of  each  layer  scales  inversely  with  B,  it  just  allows 
the  local  curvature  of  the  perpendicular  layer  and  associated  effects 
(e.g.,  gradient  and  curvature  drifts)  to  be  neglected. 

Quasi-neutrality  implies  that  the  electron  density  is  equal  to  the 
ion  density  everywhere,  i.e.,  ng  =  n^.  Because  the  electron  gyrof requency 
is  very  large  compared  to  the  other  frequencies  in  the  problem,  the 
electron  momentum  equation  perpendicular  to  the  magnetic  field  may  be 
reduced  to  the  statement  that  the  °l(j'  -  ar<-*  perfectly  magnetized  and 
move  with  the  E  x  B  velocity  within  the  ion  layer: 


The  perpendicular  ion  flow  is  not  divergence-free,  so  a  parallel 
electron  flux  or  current  is  driven  along  the  magnetic  field  between  the  ion 
layers  in  order  to  maintain  quasi-neutrality.  Ve  will  neglect  dynamic 
effects  rising  solely  from  the  parallel  dynamics,  such  as  parallel 
potential  drops,  so  that  the  parallel  electron  momentum  equation  can  be 
neglected  since  it  is  redundant  with  the  quasi-neutrality  condition.  The 
electron  and  ion  densities  are  further  affected  by  the  recombination  rate, 
but  this  has  no  effect  on  the  structure  evolution  that  we  are  interested  in 
simulating.  That  is,  recombination  acts  to  reduce  the  plasma  population, 
not  to  damp  the  development  of  structure.  As  a  result,  we  neglect 
recombination  in  our  model. 

For  each  layer  the  governing  equations  are  the  continuity  and  momentum 
equations  for  each  ion  species,  i.e.. 
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where  i  denotes  the  ion  species,  n.  the  ion  density,  v.  the  ion  velocity, 
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2^  =  q^B/ciiK  the  ion  cyclotron  frequency,  m.  the  ion  mass,  q^  the  ion 


charge,  the  ion-neutral  collision  frequency,  the  neutral  velocity,  c 


the  speed  of  light,  and  E  the  electric  field.  Temperature  effects  have 
been  neglected  in  the  ion  momentum  equations  and  all  vector  quantities  and 
gradient  operators  are  assumed  to  bp  perpend i cular  to  the  magnetic  field 


Bez .  In  the  following  ve  assume  »  Im  (  'be  i  'nosphet i •  ■  «ni ribution  to  our 
model  can  be  described  by  an  F-layer  only  in  order  to  isolate  and  study  the 


effects  of  Pedersen  conductivity  on  the  Kel vin-Helmhol tz  instability. 
Further,  we  assume  that  the  ion  gyrof requency  is  large  compared  to  both  the 
ion-neutral  collision  frequency  (F-region  approximation)  and  the 
characteristic  frequency  of  the  Kelvin-Helmhol tz  instability,  so  that  the 
ion  velocity  may  be  replaced  by  the  E  x  B  velocity  in  the  small  terms  in 
(2),  leading  to  the  following  expression  for  the  ion  velocity  when  the 
equation  is  solved  for  v . : 


v.  =  V 
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The  ion  velocity  is  therefore  equal  to  the  E  x  B  velocity  plus  small  terms 
corresponding  to  lb.:  Pedersen  drift  and  the  polarization  drift.  Note  that 
Hall  drifts  have  been  neglected  in  this  approximation  since  they  are  second 
order  in  v^/Sb.  The  ion  and  electron  densities  are  equal  by  quasi¬ 
neutrality,  and  the  ion  and  electron  velocities  are  both  equal  to  the  E  x  B 
velocity  to  lowest  order,  so  the  continuity  equation  for  each  level  has  the 
f  orm : 

I7  +  V  •  nV  =  0  where  n  =  n  =  n.  (5) 

at  -  -  e  1 


Equation  (5)  is  an  approximation  to  equation  (2),  made  by  neglecting 
Pedersen  and  polarization  drifts  relative  to  the  ExB  drift.  This  turns  out 
to  be  a  good  approximation  for  the  physical  situation  we  are  simulating, 
but  it  is  possible  to  simulate  the  full  equation  (2).  The  approximation  is 
made  here  primarily  in  the  inUac.  t  ■■!  :  ■  imp  1  i  f  7  i  ng  the  theoretical 

d i s cus  s i on . 

The  perpendicular  cuiient  d<  ir; i  t  -  ■  >  I :  j  > .  .  ...  h  1  ivei  may  be  calculated 

f 1 nm  (  \ )  and  ( 4  ) : 


where  op  -  (necv^/BS^)  is  the  Pedersen  conductivity  of  the  layer,  cM  - 

2  2  1/2  ' 
(c  /4nV^  )  is  the  inertial  capacitance,  and  ■  B/(4nnm^)  is  the  Alfven 

velocity.  In  order  to  maintain  quasi-neutrality,  the  divergence  of  the 

total  current  density  must  be  zero,  so  the  perpendicular  divergence  of  this 

current  must  be  canceled  by  a  parallel  divergence  in  the  electron  current. 

If  it  is  assumed  that  all  currents  generated  within  this  model  must  also 

close  within  this  model,  then  the  field-line  integrated  divergence  of  the 

perpendicular  current  must  be  zero,  i.e., 


0  »  /  dz  7  •  i  Jped  +  Jpol  ) 

-  »•  (tp  [<*!  VS.  ]  *  CM  [It*  i  ■  i]  s)  <7> 

where  Jped  and  Jpo^  are  the  field-line  integrated  Pedersen  and  polarization 
currents,  Ep  •  /  op  dz  is  the  field-line  integrated  Pedersen  conductivity, 
Cu  =  J  cu  dz  is  the  field-line  integrated  inertial  capacitance,  and  V  =  (J 
apvR  dz)  /  Ep  is  the  conductivity-weighted  average  neutral  velocity  along 
the  field  line.  If  we  now  define  ♦  to  be  the  electrostatic  potential, 
i.e.,  E  =  -V<fr,  the  behavior  of  the  system  is  defined  by  an  ion  continuity 
equation  at  each  level: 


|a  *  rE(V»)  -  o 


and  a  single  potential  equation  for  the  electrostatic  potential  common  to 
all  layers: 

0  -  5  •  (  [  »♦  -  !  V  S  ]  *  CM  [  §7  *  i  <8z*  ’♦>  ■  5  ]  ).  <’> 

Thus,  the  plasma  on  each  level  E  x  B  drifts  under  the  influence  of  the 
perpendicular  electric  field,  which  is  '•‘coif  determined  by  the  fact  that 
the  perpendicular  Pedersen  and  polarization  drift  currents  driven  by  the 
field  must  close  by  means  of  parallel  electron  currents  within  the  model. 


The  Pedersen  currents,  represented  by  the  first  term  in  (9),  tend  to  reduce 

any  potential  which  is  not  supported  by  the  density  gradients  and  the 

neutral  winds,  while  the  polarization  drift  currents,  represented  by  the 

second  term,  tend  to  slow  this  decay.  As  a  result,  the  balance  between  the 

effects  defines  an  inertial  relaxation  rate  \3-  E  /Cu  for  changes  in  the 

p  n 

potential. 

In  order  to  apply  (8)  and  (9),  it  is  necessary  to  specify  the  exact 
relationship  between  the  density  n  along  the  field  line  and  the  field-line 
integrated  quantities  Ep  and  C^.  In  general  this  is  accomplished  by 
separating  he  configuration  into  discrete  vertical  regions  with  similar 
characteristics,  e.g.,  a  magnet os pheric  region  and  an  ionospheric  region, 
and  treating  each  region  as  a  single  layer.  In  a  previous  paper  [Mitchell 
et  al.,  1985],  we  defined  a  raagnetospheric  layer  of  constant  density  and  an 
ionospheric  layer  with  a  density  gradient  unstable  to  the  E  x  B 

instability.  In  such  a  case,  the  magnetospheric  continuity  equation  was 
not  necessary  since  a  uniform  density  is  unchanged  by  the  incompressible 
E  x  B  flow  field.  Therefore,  it  was  necessary  to  retain  only  one 
continuity  equation,  for  the  ionospheric  density,  and  the  Pedersen 
conductivity  Ep  in  the  potential  equation  was  proportional  to  this  density. 
The  capacitance  CM,  which  is  largely  due  to  the  cumulative  inertia  of  the 
particles  along  the  long  length  of  field  line  in  the  magnetosphere,  became 
a  constant. 

In  this  paper  we  wish  to  study  the  effect  that  the  coupling  between 
magnetospheric  polarization  drift  currents  and  ionospheric  Pedersen 
currents  has  on  the  development  on  an  inertial  instability,  specifically, 
the  Kelvin-Helmholtz  instability  which  occurs  on  a  velocity  shear 
interface.  Once  again  there  is  an  iono  nheiie  layet  'ontribut ing  the 

Pedersen  conductivity  Ep  and  a  magnetospheric  layer  contributing  to  the 
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configuration.  The  velocity  shear  which  excites  the  instability  is  present 
along  the  entire  field  line  because  of  the  electric  field  mapping.  For  the 
ionospheric  layer  the  shear  configuration  requires  a  density  gradient  in 
order  to  be  an  equilibrium  in  the  presence  of  Pedersen  conductivity.  The 
source  of  the  velocity  shear  is  assumed  to  be  magnetospheric  and  naturally 
associated  with  a  density  gradient,  such  as  the  shear  at  the  plasmapause, 
so  we  will  assume  that  there  is  a  density  gradient  across  the  shear  along 
the  entire  field  line.  This  assumption  allows  us  to  view  the  entire  field 
line  as  a  single  layer  with  field-line  integrated  density  N  «  f  n  dz  and  to 
assume  that  both  the  Pedersen  conductivity  2^  and  the  inertial  capacitance 
CM  are  proportional  to  N.  Therefore  the  relaxation  rate  v  is  a  constant 
for  this  model,  and  the  effect  of  the  Pedersen  currents  on  the  development 
of  the  Kelvin-Helmholtz  instability  depends  on  the  relationship  of  the 
value  of  this  rate  to  that  of  the  growth  rate  of  the  instability.  If  v  is 
much  smaller  than  the  growth  rate,  the  model  is  purely  inertial,  and  if  the 
two  rates  are  comparable,  the  collisional  effects  are  significant.  The 
explicit  form  of  the  continuity  and  potential  equations  for  this  model  now 


take  the  form 


[  It  *  I  <v  ”♦>  ■  7  ] N 


0  =  v 


( ;  m  [  »♦  -  “  v  'z  ]  *  N  L  It  *  I  7+) 


;  »♦)_ 


where  explicit  use  of  the  incompressibility  of  the  E  x  B  flow  has  been  used 
in  (10). 

The  basic  equililibrium  configuration  for  the  Kelvin-Helmholtz 


instability  consists  of  a  •eloci 


’  I  7  )  o 


upported  by  an 


electric  field  E  =  E^(y)  such  that  V^fy)  =  cEy(y)/B. 


*-■  V*  N  -  , 

na'vvv.n-; 


\ vs’ 


There  are  no 


polarization  currents  in  this  configuration  because  the  gradient  of  E  is 


perpendicular  to  V.  If  there  are  no  ion-neutral  collisions,  this 


configuration  will  be  an  equilibrium  for  any  density  profile  for  which 


3N/3x  *  0.  In  the  presence  of  ion-neutral  collisions,  the  Pedersen 


conductivity  is  non-zero  and  there  are  Pedersen  currents  in  the  y- 


direction,  so  that  this  configuration  is  an  equilibrium  only  if  the 


divergence  of  the  Pedersen  current  is  zero  everywhere,  i.e.,  if 


3N(y)Ey(y)/3y  =•  -  V*(N(y)Vn  x  B/c).  If  we  make  the  simplifying  assumption 


that  the  neutral  wind  is  uniform,  the  equilibrium  condition  becomes 


3(N(y)(Ey(y)  -  VnxB/c))/3y  -  0.  This  condition  reflects  the  fact  that  the 


Pedersen  currents  act  to  short  out  the  electric  field  associated  with  the 


velocity  shear  on  the  time  scale  \T  unless  the  field  is  supported  by  the 


current  due  the  neutral  collisions  with  the  ions.  Although  the  equilibrium 


allows  an  arbitrary  neutral  wind  component  parallel  to  the  density 


gradient,  we  will  assume  that  V  *  0,  since  a  parallel  neutral  component 


introduces  E  x  B  instability  effects  into  the  configuration.  An 


investigation  of  the  E  x  B/Kelvin-Helmholtz  instability  coupling  for  an 


arbitrary  neutral  wind  direction  will  await  a  later  work. 


In  our  calculations  we  explicitly  assume  that  the  equilibrium  velocity 


shear  has  the  form  V  (y)  =  -Vn  tanh  (y/L)  so  that  the  shear  in  V  is 


localized  within  a  scale  length  L  of  the  origin  in  y,  and  is  roughly 


constant  outside  of  this  region.  This  form  implies  that  the  equilibrium 


electric  field  and  density  have  the  forms  E^(y)  =  -(BV^/c)  tanh  (y/L)  and 


N(y)  =  NqVq/{  -  Vq  tanh  (y/L)  -  ] .  Note  here  that  the  magnitude  of  Vn 


must  be  greater  than  that  of  Vn  or  rhete  will  be  a  singularity  in  N(y). 


The  implication  is  that  there  must  b«  i  o^n-'err>  neutral  wind  in  the  frame 


in  which  the  velocity  shear  is  symm<  1 i ' 


.plicit  -h"ice  of  the  ratio 


IWVq  determines  the  magnitude  of  the  density  jump  across  the  velocity 


,  t1*  **»**»» 
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shear,  i.e.,  the  value  of  N(  +  °°)/N{ -°°) .  In  the  cases  we  will  show  here  this 


ratio  is  taken  to  be  -2  so  that  N(+®)/N(-“)  =  3.  Another  point  is  that  the 


anti-symmetric  nature  of  the  Kelvin-Helmholtz  equilibrium  in  the  y- 


direction  has  been  broken  by  the  required  form  of  the  density  gradient. 


The  geometry  of  the  equilibrium  we  have  chosen  is  shown  in  Fig.  2.  This 


particular  equilibrium  is  chosen  since  it  simulates,  approximately,  the 


plasma  configurations  along  the  plasma  sheet  outer  and  inner  boundaries  and 


its  footprints  in  the  polar  ionosphere  in  and  near  auroral  and  polar  cap 


It  is  useful  to  note  here  that  the  equilibrium  geometry  that  we  have 


defined  occurs  naturally  in  the  non-linear  development  of  the  E  x  B 


instability.  In  that  instability,  a  neutral  wind  blowing  down  a  density 


gradient  causes  the  gradient  to  striate  in  the  direction  perpendicular  to 


the  neutral  wind  and  'fingers'  to  grow  out  of  the  gradient.  The  edges  of 


these  fingers  have  the  equilibrium  geometry  defined  here,  i.e.,  a  velocity 


shear  and  associated  density  gradient  perpendicular  to  a  neutral  wind.  In 


our  previous  paper  [Mitchell  et  al.,  1985],  it  was  shown  that  the  inclusion 


of  inertial  effects  in  this  configuration  caused  the  growth  of  modes  with 


wave  vectors  parallel  to  the  neutral  wind  direction,  a  result  due  possibly 


to  Kelvin-Helmoltz-like  effects  in  the  finger  geometry. 


III.  LINEAR  THEORY 


In  order  to  study  the  linear  theory  of  the  growth  of  the  Kelvin- 


Helmoltz  instability  within  the  equilibrium  we  have  defined,  the  linear 


perturbation  is  assumed  to  have  wave  number  k  in  the  x-direetior.  and 


frequency  w,  so  that  the  quantifier  rl  "i-i  1  ""v  be  written: 


N  =  N(y)  +  6N(y)  el(kx"wt)  ;  <f>  -  *(y)  +  8«Ky)  el(kx_Wt)  .  (12) 


•  *  v  v.  f.  .*  • 


By  substituting  these  expressions  into  (10)  and  (11),  we  obtain 


(» -  kv  ®  f  |  *♦ 


(»  *  i5  -  KVX)  (f^M)  -  Hk2S+) 

-  -  k“  fc("  b  vx)  -  l5l  *(»„  -,vx)  57(f)  •  <14> 

Equation  (13)  can  be  substituted  into  (14),  yielding  the  following  mode 
equation  for  the  Kelvin-Helmholtz  instability 

hi'hi^rj]  -  kM^v  ) 


where 

(w  -  kV  )  («  +  iv  -  kV  ) 

*  (V  -  V  ) 

v  n  x' 

For  the  region  outside  the  shear  layer,  Vx  and  A  are  constant  so  that  (15) 


reduces  to 


^~2  5*  =  k26*  , 


which  implies  that  the  perturbation  potential  has  a  decay  length  of  1/k  in 
the  shear  direction. 

We  have  solved  (15)  analytically  (see  Appendix  A)  in  the  sharp  boundary 
limit  and  numerically  for  a  smooth  boundary  profile.  For  the  smooth 
profile  we  have  taken  =  -V^  tanh  (y/L)  and  the  results  are  plotted  in 
Fig.  3  for  the  cases  \>  =  \>/(V  /L'  ' ' .  M  1.  n.:,4.  o.5.  The  'classic' 
Kelvin-Helmholtz  case  is  represented  by  the  upper  curve  in  this  figure,  and 
is  obtained  by  setting  v  =  0  and  assuming  Vn  >>  in  (15).  The  latter 


assumption  reduces  the  background  density  gradient  to  zero.  In  this  case 


the  growth  rate  y  (where  «  -  +  iy)  has  a  maximum  value  of  y  «  0.19  Vq/L 
at  kL  *  0.44.  For  k  less  than  this  value,  the  growth  rate  is  reduced 
because  the  width  of  the  shear  region  L  is  small  compared  to  the  decay 
length  1/k  of  the  potential  in  the  shear  direction,  and,  for  k  greater  than 
this  value,  the  growth  rate  is  reduced  because  the  effective  shear  within 
the  length  1/k  is  small.  This  case  is  symmetric  in  the  shear  direction  so 
the  phase  velocity  wr/k  of  the  instability,  is  zero.  The  addition  of  a  3:1 
density  gradient  across  the  velocity  shear  (V^Vq  m  2)  reduces  the  maximum 
growth  rate  by  about  102,  and  produces  an  asymmetry  in  the  shear  direction, 
resulting  in  a  real  phase  velocity  u^/k  *  -.13  Vq.  A  10:1  density  gradient 
(Vo/Vq  =•  1.22)  reduces  the  maximum  growth  rate  by  50X  [ Satyanarayana  et 
al.,  1987b]. 

The  inclusion  of  Pedersen  conductivity  further  reduces  the  growth  rate 
of  the  instability.  From  (15)  it  can  be  seen  that  the  growth  rate  of  the 
instability  depends  solely  on  the  behavior  of  the  function  A.  This 
function  may  also  be  written 

A  .  [(»  .  |  5  -  kVx)2  -  \  V2  ]  <Vn-  Vx)-‘  (17) 

If  v  is  small  compared  to  kVx,  v  is  negligible  in  this  equation,  and  the 
growth  rate  is  reduced  by  a  factor  of  v/2  from  the  collisionless  case. 
This  can  be  seen  in  the  plots  of  growth  rate  for  the  collisional  cases  v  s 
\J/(Vo/L)  =  0.10,  0.34,  and  0.50  in  Fig.  3.  The  growth  rate  is  reduced  by  a 
significant  amount  when  v  becomes  comparable  to  the  collisionless  growth 
rate  (about  .2  Vq/L),  and  the  instability  is  effectively  damped  for 
v  >  0.5  Vq/L.  As  the  density  gradient  is  increased,  the  growth  rate  for 
all  of  the  cases  is  decreased.  Ve  find  good  agreement  between  the  growth 


rate  as  derived  from  the  sharp  boundary  analysis  in  the  small  wavenumber 
limit,  as  outlined  in  Appendix  A,  and  that  computed  numerically  from  Eq. 


IV.  NUMERICAL  SIMULATIONS 

In  our  numerical  simulation  of  the  evolution  of  the  Kelvin-Helmholtz 
configuration  that  we  have  described,  the  calculation  is  performed  on  a  64 
x  100  cell  grid  (x,y)  which  is  periodic  in  the  x-direction  and  has 
Dirichlet  boundary  conditions  in  the  y-direction  (potential  set  equal  to  a 
constant  at  the  boundaries).  In  the  x-direction,  he  cell  size  is  uniform 
and  contains  or.e  wavelength  of  the  unstable  perturbation  for  the  Kelvin- 
Helmholtz  instability  within  the  64  cells.  The  y-direction  cell  size  is 
uniform  tor  the  80  cells  covering  the  interior  region  from  y  -=  -2L  to  +2L, 
with  the  10  boundary  cells  on  either  end  being  a  stretched  mesh  to  y  =  -13L 
and  el3L.  The  aspect  of  the  cells  in  the  uniform  mesh  region  is  4:1  for 
the  wavelength  used.  The  numerical  method  used  to  simulate  (JO)  and  (11) 
is  described  in  Zalesak  et  al.  [1982].  The  continuity  equation  is  solved 
numerically  using  the  multidimensional  flux-corrected  transport  techniques 
of  Zalesak  [1979],  while  the  potential  equation  is  solved  with  the 
incomplete  Cholesky  conjugate  gradient  algorithm  of  Hain  [1980]. 

We  have  performed  three  Kelvin-Helmholtz  simulations  with  this  model, 
all  of  which  have  a  3:1  density  jump  and  a  wavenumber  of  kQ  =  0.44/L.  The 
three  cases  have  scaled  inertial  relaxation  tates  of  v  0.00.  0.]n,  arid 
0.34,  and  are  represented  by  the  dots  on  the  linear  growth  rate  curves  of 
Fig.  1.  These  three  cases  have  1  itv'-i*  ft  rates  ><f  v.  <v  .  !.)  ■  i .  j  ' , 


o.l2.  and  O.Oi,  respec t  i  v<  ly .  vhicit 


tin  maximum  growth 


rates  for  their  respective  values  of  v.  The  simulations  are  initialized 
with  the  appropriate  eigenfunction  for  this  mode,  determined  by  a  numerical 


e  i  gent  tine  t  ion  code,  so  that  the  !  ami:  hehav  i< 


>t  the  t vo -d i mens i ona 1 


elect! nstai ic  code  may  be  verified. 


!i«  •• :  gent ;  1 1 .  *  t  ions  tor  these  three 


cases  are  plotted  in  Fig.  *4.  The  ampli'ude  and  phase  of  the  Fouriet 
components  6*  are  shown  in  Fig.  aa  and  ah.  respectively,  and  the  full  2  - 


dimensional  eigenfunctions  are  contoured  in  Fi; 


For  each  case  the 


eigenfunction  has  a  phase  shift  within 


the  siieai  region  about  tne  origin. 


w  w  m  •xrvmmmm wwro 


with  the  total  phase  shift  increasing  from  about  1  radian  for  v  .  0.00,  to 
about  1.6  radians  for  v  -  0.34. 

The  behavior  of  the  densities  and  potentials  for  the  three  simulations 
are  displayed  in  Figs.  5-7.  In  the  plot  of  the  \>  *  0.00  simulation,  the 
density  gradient  serves  primarily  as  a  tracer  for  a  standard  Kelvin- 
Helmholtz  simulation.  The  initial  shear  boundary  perturbation  grows  in 
amplitude  and  width  until  the  width  is  on  the  order  of  the  shear  scale 
size,  at  which  time  the  boundary  layer  forms  the  characteristic  vortices 
and  spirals  of  the  Kelvin-H  Imholtz  instability  (shown  in  the  last  panel  of 
Fig.  5a).  In  our  simulation,  numerical  effects  begin  to  diffuse  small 
scale  structure  when  the  width  of  the  boundary  spiral  approaches  the  cell 
size.  This  occurs  at  roughly  time  8t  in  our  simulation. 

The  v  -  0.10  simulation  is  shown  in  Figs.  6a  and  6b.  The  linear  phase 
of  this  simulation  is  similar  to  the  \>  «  0.00  case.  When  the  vortices 
begin  to  form,  however,  the  Pedersen  conductivity  acts  to  prevent  curling 
of  the  density  contours,  and  the  contours  trail  out  from  the  perturbation 
as  the  perturbation  convects  to  the  left.  Since  the  structure  is  now 
elongated  in  the  x-direction,  it  becomes  unstable  to  the  Kelvin-Helmhol tz 
instability  on  this  new  interface.  Smaller  -^ale  instabilities  due  to  this 
effect  are  evident  in  the  last  panel  tic  a.  In  addition,  the  leftward 
convecting  perturbation  structure  overtakes  the  trailing  density  structure, 

P(k*>  ■  E  <l£k  ,k  lz> 

k  *  y 
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p<k  >  ■  E  <|fk  k  l2> 

k  *  y 

x 

We  now  give  selected  power  spectra  for  the  cases  v  =  0.00,  0.10,  and  0.34. 

The  typical  temporal  evolution  of  the  potential  Fourier  modes  for 
several  values  of  k^  is  shown  in  Fig.  8  for  the  v  =  0.00  case.  The 
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amplitude  of  the  dominant  growing  mode  k^  =  k^  increases  exponentially 
until  nonlinear  effects  become  important  and  the  mode  reaches  a  saturation 
level.  Agreement  between  the  theoretical  value  of  the  linear  growth  rate 
and  measured  value  is  excellent.  In  all  three  cases,  \>  »  0.00,  0.10,  0.34, 

our  simulations  match  the  predicted  linear  growth  rate  for  2  to  4  growth 

times  after  initialization.  The  higher  harmonics  (2kg,  3k^,  etc.)  increase 
in  amplitude  at  a  faster  rate  until  saturation  since  they  are  driven  by  the 

dominant  mode.  At  saturation  the  modes  oscillate  in  amplitude  as  they 

exchange  energy  with  each  other. 

Figures  9  and  10  give  the  time  and  transverse  averaged  power  spectra  of 

the  potential  5$  =  ♦  -  ♦q  and  density  &N  =  N  -  N^,  respectively,  for  the 

case  v  -  0.00.  The  wavenumbers  k  and  k  are  in  units  of  k~  =  2n/L  = 

x  y  Ox  x 

0.5/L  and  kn  *  2n/L  =  1.6/L,  respectively,  where  L  is  the  length  of  the 
uy  y  x 

system  in  the  x  direction,  Ly  is  the  length  of  system  in  the  uniform  mesh 

region  in  the  y-direction,  and  L  is  the  half-width  of  the  shear  layer.  We 

find  P(k  )  «  k  ~5*4  4  °'4  and  P(k  )  «  k *5'3  4  °'3  for  |&<f>|2  and  P(k  )  - 

xx  y  y  x 

k  -2‘2  4  0,2  and  P  ( k  )  k  ~2'3  4  °'3  for  1 5N  | 2 .  The  spectra  P(k  )  for  the 
x  y  y  x 

potential  &<fr  is  in  approximate  agreement  with  that  computed  in  the  work  of 
Miura  and  Sato  (1978)  which  showed  |5<Kk...v  t  kw  .  Our  measured 

and  multiple  overlapping  structures  are  created  with  smaller  y-direction 
scales  rl'.an  in  the  previous  case.  The  resulting  non-linear  structure  of 
the  instability  is  significantly  different  from  the  v  =  0.00  simulation, 
even  rhough  the  linear  growth  rates  differ  only  by  about  30  percent. 

The  nonlinear  behavior  of  the  v  =  0.34  case  is  again  different  from  the 
two  pi  oreed  ing  cases.  In  this  Simula  f  t  on  .  d<*n«:i  ty  contours  are  inhibited 
from  curling  very  early  in  the  development  ol  the  instability,  and  so  the 
y-direction  extent  of  the  perturbation  is  reduced.  Density  is  swept  off 
the  high  density  peak  and  opposite  to  the  direction  of  the  phase  velocity, 
resulting  very  quickly  in  a  smaller  scale  sttucture  parallel  to  the  shear 


direction  than  in  either  of  the  other 


■;  i  mu  i  .i  ’  i  ons  .  The  density  plots  in 


Fig.  7a  clearly  show  the  increased  small-scale  structure  and  reduced  y- 
extent  in  the  nonlinear  development  of  the  instability.  In  Fig.  7b,  the 
potential  plots  exhibit  only  small  deviations  from  equilibrium  relative  to 
the  behavior  in  the  other  tvo  simulations. 


V.  POWER  SPECTRA 


We  have  computed  the  spatial  power  spectra  of  the  density  and  potential 
fluctuations  associated  with  the  evolution  of  the  Kelvin-Helmholtz 
instability  both  with  and  without  Pedersen  coupling. 

We  consider  a  two-dimensional  spatial  Fourier  representation 


f(x,y,t) 


£  fk  ,k (t)  exp[i(kxx  +  y)] 

l  u  x’  y  3 


with  f  a  N,<fr.  Here,  f^  ^  are  the  Fourier  coefficients  defined  by 

x’  y 

fk  k  =  t2*2')1  f  dxdY  f(x-v)  “"pf-  *  K,x  +  ktryl] 
x9  y  **  “  ** ' 

The  transverse  averaged  power  spectra  P(k  )  and  P(k  )  are  defined  by 

^  2 

P(kx)  for  is  shallower  since  we  are  essentially  integrating  |5$|  over 

ky.  In  addition,  Miura  and  Sato's  spectral  results  are  not  time  averaged 

and  are  obtained  using  different  initial  conditions. 

In  Figs.  11  and  12  we  give  the  power  spectra  of  the  perturbation 

potential  and  density  &N,  respectively,  for  the  case  v  =  0.10.  We  have 

time-averaged  the  spectra  computed  in  the  nonlinear  regime  over  a  time 

interval  of  approximately  five  linear  growth  times  after  saturation  has 

2 

occurred.  The  two-dimensional  power  spectra  |f^  ^  |  ,  with  f  «  N,$,  are 

x’  y 

computed  at  twenty  different,  successive  times  in  this  interval  and  then 

averaged.  This  averaging  minimizes  numerically  induced  fluctuations  and 

gives  a  more  accurate  representation  of  the  turbulence  spectrum  in  the 

2  4  7-1-02 

saturated,  nonlinear  regime.  We  obtain,  for  |S$|  ,  P(kx)  a  kx 


and  P(k  )  «  k  ~5'3  4  0-2  while  for  |&N|2,  P(k)  «  k v'1,9  1  0,2  and  P(k  )  « 
y  y  xx  y 

kyl  **  ±  0*2_  ye  note  that  these  spectra  are  shallower  compared  to  the 
previous  case  with  v  =  0.00. 

Figures  13  and  14  display  the  power  spectra  of  6$  and  SN,  respectively, 

for  v  =»  0.34.  For  S$,  we  find  P(kx)  *  *  ®‘2  and  P(k^)  « 

ky_5‘3  *  °’2-  For  *»  we  obtain  P(kx)  «  kx"2,1  ±  °‘2  and  P(k  )  «  k  ~1,7  4 
0  3 

.  These  spectral  slopes  are  also  shallower  than  the  v  *  0  case  as  shown 
in  Figs.  9  and  11. 

The  spectral  indices  of  the  power  spectra  of  the  total  electrostatic 
potential  <fr  and  density  N  are  listed  in  Table  I  and  display  similar  values. 

In  general,  we  find  shallower  power  spectra  for  the  v  =  0.10,  0.34 
cases  as  opposed  to  the  \>  =  0-00  case.  Ue  feel  that  this  effect  is  caused, 
in  part,  by  the  excitation  of  shorter  wavelength  (higher  k)  potential 
fluctuations  by  secondary  Kelvin-Hel  mh''l  *  in-tabilities.  as  noted  in  Sec. 
IV,  in  the  nonlinear  regime  thus,  increasing  the  power  at  higher  k,  and,  as 

a  result,  shallowing  the  spectra.  We  give  a  simple  argument  for  this 
shallowing  of  the  spectra  based  on  the  lowest  oraer  nonlinear  evolution  of 
the  Kelvin-Helmholtz  instability  in  Appendix  B.  Also,  in  addition  to  the 
enhanced  power  at  short  wavelengths,  we  note  that  there  is  a  strong 
reduction  of  power  at  long  wavelengths,  i.e.,  small  k,  for  the  collisional 
cases.  For  example,  there  is  an  order  of  magnitude  reduction  in  the  power 
at  k^2  =  4  for  v  =  0.34  (Fig.  16b)  from  v  =  0.00  (Fig.  12b).  This  is 

because  the  large-scale  vortices  (or  structures)  which  form  in  the  inertial 
regime  (\>  =  0.00)  and  cause  the  enhanced  power  at  long  wavelengths  are 
suppressed  in  the  collisional  regime  ( \>  -  0.3M  which  therefore  leads  to  a 

decrease  in  the  power  at  long  wavelengths. 

VI.  .SUMMARY  AND  CONCLUSIONS 


In  this  paper  we  have  studied,  using 
simulation  techniques,  the  linear  and 


both  analytical  and  numerical 
nonlinear  evolution  of  the 


V* 


electrostatic  Kelvin-Helmholtz  instability,  resulting  from  velocity-sheared 

plasma  flows  perpendicular  to  a  magnetic  field,  including  Pedersen 

conductivity  effects.  We  find  that  the  Kelvin-Helmholtz  instability 

develops,  in  the  nonlinear  regime,  in  a  distinctly  different  manner  when 

Pedersen  conductivity  effects  are  included  than  the  case  where  there  is  no 

Pedersen  conductivity  (i.e.,  no  ion  collisions).  We  find  that  Pedersen 

coupling  effects,  in  conjunction  with  a  neutral  wind  and  plasma  density 

gradient,  (1)  lead  to  smaller  growth  rates  (increased  time  scales),  in  the 

linear  phase  of  the  Kelvin-Helmholtz  instability,  (2)  inhibit  Kelvin- 

Helmholtz  vortex  formation,  (3)  lead  to  nonlinear  structures  which  can  be 

described  as  "breaking  waves",  and  (4)  generate,  in  the  nonlinear  regime, 

small  scale  structure  and  turbulence  by  means  of  secondary  Kelvin-Helmholtz 

instabilities  growing  on  primary  Kelvin-Helmholtz  driven  waves.  Ue  have 

also  computed  the  spatial  power  spectra  of  the  electrostatic  potential  and 

density  fluctuations  associated  with  the  Kelvin-Helmholtz  instability  and 

find  differing  spectra  with  and  without  Pedersen  conductivity  effects.  For 

the  standard  Kelvin-Helmholtz  velocity  sheared  configuration  V  =  Vx(y)  ex, 

but  with  a  density  gradient  n^(y)  perpendicular  to  the  velocity  flow,  we 

find  P(k  )  «  k  ±  and  P(k  )  k  for  the  electrostatic 

xx  y  y 

potential  fluctuations,  and  P <k^ )  *  k^-^ ‘ ^  t  and  P(k^)  *  ky-^'^  ± 

for  the  density  fluctuations,  where  P(k  )  and  P(k  )  are  the  time-averaged, 

x  y 

spatial  power  spectra  averaged  over  the  y-  and  x-directions,  respectively. 

For  the  Kelvin-Helmholtz  instability  with  a  density  gradient  but  including 

ion-neutral  collisional  effects  (Pedersen  coupling)  we  find,  in  general, 

shallower  spectra  with  P(k  )  *  k  ~(^*^-5)  ancj  p(k  )  <*  ^  -(5-5.5)  j  the 
r  '  x  x  y  y 

potential  fluctuations  with  P(k^)  *  k^  ‘•■5)  ancj  P(k^)  11  for 

the  density  fluctuations. 

Recently,  Basu  et  al.  [1987],  using  Dynamics  Explorer  satellite  data, 
have  measured  the  simultaneous  density  and  electric  field  fluctuation 
spectra  associated  with  velocity  shears  in  both  the  southern  and  northern 


auroral  oval.  In  this  region  they  observe  strong  shear  when  E„  >  E  and 

H  ~  p 

moderate  shear  when  Ey  <  E^  where  Ey(Ep)  is  the  approximate  integrated  Hall 

(Pedersen)  conductivity.  In  the  moderate  velocity  shear  region  they  find 

shears  V/L  =  1  sec-1  (where  V  -  1  km/sec  and  L  -  1  km),  and  Pedersen 

conductivities  I  “1-5  mhos.  Taking  a  typical  C  =10  farad  where  C  is 
p  m  m 

defined  in  Section  2  we  find  v  -  ( Ep/Cm)/(VQ/L)  =  0.1  -  0.5  which  includes 

the  values  selected  for  v  in  Section  3  and  4.  In  the  velocity  shear  region 

2 

Basu  et  al.  [1987]  compute  the  transverse  averaged  power  spectra  |E(kg)  |  « 

kg  P,  p  =  2.8  -  3.6  and  |n(kg) *  kg  p  =  1.7  -  2  where  E(n)  denotes  the 

electric  field  (density).  Here  s  denotes  rhe  direction  of  the  satellite  as 

it  passes  through  the  velocity  sheared  rer.  i-n.  Taking  E(k)  <*  k<j>(k)  this 

gives  for  the  potential  fluctuations  |$(k  )|^  *  k^-15,  p  =  4.8  -  5.6.  It  is 

difficult  to  relate  k  to  the  coordinate  system  used  in  the  simulations  of 

— s 

Sec.  4  since  the  initial  conditions  for  the  assumed  Kelvin-Helmholtz 

fluctuations  in  the  ionosphere,  as  observed  by  satellites,  are  not  known. 

The  observations  of  Basu  et  al.  [1986]  show  velocity  shear  gradients  along 

the  satellite  track  suggesting  that  k  =  k  (direction  of  velocity  gradient 

s  y 

at  t  *  0  in  the  simulations  of  Sec.  4).  Taking  k  »  k  the  observed 

s  y 

spectral  indices  for  the  electric  field  (potential)  and  density  are  in 
reasonable  agreement  with  those  computed  from  the  simulations  of  Section  4 
where  we  found,  for  v  =  0.10,  0.34,  P(k  )  «  k  5.*+)  £  the 

y  y 

potential  and  P(k^)  «  ky  2.0)  for  t^e  dens  ity.  Our  results  are  for 

k  and  k  in  the  range  0.5  -  5  km*^  and  1.6  -  16  km  *,  respectively,  or  for 
x  y 

scale  sizes  in  the  range  12.4  km  -  1.2  km  and  3.9  -  0.39  km,  respectively. 

Exact  agreement  can  not  be  expected  since  ve  have  not  included  Hall 
conductivity,  parallel  potential  drops,  and  three-dimensional  effects  in 
our  simulations. 


On  larger  spatial  scales  our  simulation  results  approximate  plasma 
conditions  along  the  magnetospher ic  tail  plasma  sheet  outer  and  inner 
boundaries,  and  its  footprints  in  the  polar  ionosphere,  the  auroral  oval. 


Figure  15  is  a  sketch  of  these  regions.  At  the  poleward  boundary  of  the 
plasma  sheet  -  auroral  oval  velocity  shear  is  developed  between  the 
tailward  convection  of  plasma  across  the  polar  cap  and  the  sunward 
convection  at  lower  latitudes  in  the  auroral  oval  and  the  inner  plasma 
sheet.  There  is  also  a  strong  density  gradient  across  the  shear  boundary 
from  the  very  low  densities  of  the  ionospheric  polar  cap  and  magnetic  tail 
lobe  to  the  increased  density  in  the  ionospheric  auroral  plasma  and  the 
plasma  sheet.  The  section  of  the  boundary  A  -  A'  shown  in  Fig.  15  (panels 
a  and  b)  is  similar  to  the  shear  layer  model  in  this  simulation. 

The  section  of  the  boundary  B  -  B'  shown  in  Fig.  15  (panels  a  and  c) 

also  exhibits  a  velocity  shear  layer.  Here  the  shear  is  developed  between 

the  sunward  convection  in  the  auroral  oval  and  plasma  sheet  and  the  slowly 

convecting  plasma  at  lower  latitudes  in  the  ionosphere  and  closer  to  the 

earth  in  the  magnetotail.  Again  in  these  regions  a  density  gradient  exists 

across  the  shear  layer  from  low  latitude  plasma  toward  higher  latitude 

plasma.  In  order  to  assess  the  efficiency  of  the  Kelvin-Helmholtz 

instability  for  the  plasma  sheet  -  auroral  region  boundary,  it  is  necessary 

to  estimate  the  quantities  Cm  and  in  these  regions.  For  a  moderately 

disturbed  nightside  auroral  oval  the  conductivity,  E^,  is  typically  5-10 

mhos.  For  the  outer  plasma  sheet  we  estimate  C  to  be  about  600  farads, 

m 

_3 

taking  the  density  to  be  1  cm  ,  the  average  field  strength  to  be  5  nT,  and 
a  typical  integration  length  of  10  km.  For  the  inner  plasma  sheet  is 

_3 

considerably  smaller  and  about  10  farads,  taking  the  density  as  1  cm  ,  the 

4 

field  strength  as  30  nT,  and  again  a  typical  integration  length  of  10  km. 
For  shear  frequencies  Vq/L  -  0.1  -  1.0  sec  ^  the  quantities  estimated  above 
would  indicate  that  the  high  latitude  auroral  zone  -  plasma  sheet  shear 
boundary  should  be  unstable  to  the  Kelvin-Helmholtz  shear  instability  with 
only  moderate  modifications  by  the  ionospheric  conductivity  because  v  < 
0.1.  On  the  other  hand  the  low  latitude  auroral  -  inner  plasma  sheet  shear 
boundary  should  be  strongly  affected  by  the  ionospheric  conductivity 


because  v  >  0.5,  and  it  may  be  only  marginally  unstable  depending  on  the 
precise  balance  between  plasma  inertia  and  conductivity.  Figure  16,  a 
sketch  of  the  auroral  morphology  from  Akasofu  [1976],  shows  this  general 
behavior  with  a  very  highly  irregular  and  structured  poleward  boundary  and 
a  much  less  disturbed  boundary  at  lover  latitude. 

Finally,  we  wish  to  point  out  that  our  simulation  model  for  simplicity 
has  not  included  Hall  conductivity,  parallel  potential  drops  and  associated 
imperfect  magnetospheric-ionospheric  coupling,  variations  along  the 
geomagnetic  field,  and  recombination  terms.  As  a  consequence,  our  results 
are  not  strictly  applicable  to  the  auroral  zone  plasma.  It  was  our 
intention  to  isolate  and  study  only  the  effects  of  Pedersen  conductivity  on 
the  Kelvin-Helmholtz  instability  as  it  may  occur  in  the  high  latitude 
ionosphere  and  inner  magnetosphere.  It  is  clear  that  effects  such  as  Hall 
conductivity  and  parallel  electric  fields  are  important  in  the  auroral 
regions,  and  we  intend  to  study  their  effects  on  the  evolution  of  Kelvin- 
Helmholtz  instability  in  space  plasmas  with  more  realistic  and  self- 
consistent  magnetosphere-ionosphere  coupling  in  future  studies. 
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compare  the  analytical  results  with  the  exact  numerical  solutions  of  (Al); 
the  numerical  solutions  pertain  to  flows  with  smooth  profiles  (hyperbolic 
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From  (A14)  we  can  see  that  the  analytical  growth  rate  is 
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where  L(6<f>)  represents  terms  linear  in  6<J>.  The  expression  for  V  represents 

the  nonlinear  interactions  between  modes  k'  and  k"  such  that  k'  +  k"  =  k. 

The  first  two  terms  in  the  expression  for  V  give  the  approximate  wave 

coupling  proportional  to  k  x  k"  •  k'  when  \>  =  0.00.  This  term  vanishes  if 

k'  =  k^"  x  and  k'  =  k  '  x.  As  a  result,  in  this  lowest  order  approximation 

for  nonlinear  coupling,  the  mode-mode  interactions  along  the  k^  direction 

is  enhanced  when  v  >  0  since  its  coupling  coefficient  *  k"  •  k'  *  0  i f  k  " 

-  -  -x 

3 

-  k^"  x  and  k'  =  k  '  x.  Furthermore,  this  coupling  ( v  >  0)  has  V  k  and, 
thus,  favors  higher  wavenumbers.  This  is  consistent  with  the  spectral 
results  computed  in  Sections  4  and  5  which  indicated  small  scale  (high  k) 
structure  generation  when  v  >  0  in  the  nonlinear  regime. 


A  diagram  of  the  Kelvin-Helmholtz  shear  configuration  in  the  presence 
of  magnetosphere-ionosphere  coupling.  The  dashed  line  represents  the 
unstable  shear  interface,  and  the  instability-driven  currents  are 
indicated:  the  perpendicular  ion  polarization  and  Pedersen  currents  in 
the  magnetosphere  and  ionosphere  respectively  and  the  electron  closing 
current  parallel  to  the  magnetic  field. 


The  linear  growth  rates  of  the  Kelvin-Helmhol tz  instability  versus 
wave  number  k  for  the  'classic'  (no  density  gradient)  case  and  for 
four  collisional  cases  (v  =  0.0,  0.10,  0.34,  0.50)  with  a  density  jump 
of  3:1.  The  bullets  mark  the  three  simulation  cases. 


Fig.  4  --  Eigenfunction  plots  for  the  three  simulation  cases,  a)  The  amplitudes 
of  the  Fourier  components  5<f>  versus  y.  b)  The  phases  of  the  Fourier 
components  S<J>  versus  y-  c)  Contour  plots  of  the  complete  2-dimensiona] 
eigenfunctions.  Dashed  lines  represent  negative  contours. 
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Fig.  4 (Cont’d) —  Eigenfunction  plots  for  the  three  simulation  cases,  a)  The  amplitudes 
of  the  Fourier  components  6$  versus  y.  b)  The  phases  of  the  Fourier 
components  6<f>  versus  y.  c)  Contour  plots  of  the  complete  2-dimensional 
eigenfunctions.  Dashed  lines  represent  negative  contours. 


eigenfunctions.  Dashed  lines  represent  negative  contours. 


6.88t  where  t  =  \>  =  6.0  L/V  .  a  >  Density,  b)  Potential  cont 


POTENTIAL =  nu 
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7t 

Fiji.  X  Plot  of  amplitudes  of  tiie  four  most  significant  Fourier  components  of 
the  potential  l'f’(cx(t)|  vs<  for  r*10  v  “  0.00  case. 


4X 


Plot  of  the  time  and  transverse  averaged  power  spectra  (a)  P(k  )  and 
(b)  P(k  )  for  6N  for  v  =  0.00.  The  units  aie  described  in  Fig.  9. 
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Fig  1.1  —  Plot  of  the  ti 


(b)  P(ky)  for 


Sketch  of  regions  of  strong  velocity  shear  in  high  latitude  space 
plasma,  showing  the  (a)  polar  ionosphere  viewed  from  above,  (b) 
magnetotail  as  viewed  toward  the  sun,  and  (c)  equatorial  plane  and 
plasma  sheet  of  the  magnetosphere. 
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C  0  M  M  A  N  l>  E  h 

U.S.  ARMY  MISSILE  INTELLIGENCE  AGENCY 
REDSTONE  ARSENAL,  AL  35809 
JiC.  Y  ATTN  JIM  GAMBLE 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  DC  20375 


0  1  C  Y 

ATTN 

CODE 

44  700 

S . L .  Osuakow  , 

2  6  C 

YS  IF 

UNCLASS 

(  01C 

Y  IF 

CLASS ) 

ATTN 

CODE 

44  7  80 

J.D.  HUB  A,  CO 

r  v  c; 

IE  UNCLASS, 

0 1  C  Y  IF  CLAP.: 

0  1CY 

ATTN 

CODE 

44  7  0  1 

I.  V ITKOV  ITSK 

01  C  Y 

ATTN 

CODE 

7500 

0  1  C  Y 

ATTN 

CODE 

7  5  50 

01  CY 

ATTN 

CODE 

7580 

0  1C  Y 

ATTN 

CODE 

7551 

01  CY 

ATTN 

CODE 

7555 

0  1CY 

ATTN 

CODE 

14  7  30 

E.  MCLEAN 

01  C  Y 

ATTN 

CODE 

'4  7  5  2 

0  1  C  Y 

ATTN 

CODE 

14  7  30 

B.  R I P  I  N 

20C  Y 

ATTN 

CODE 

2  62  8 

COMMAN OK 

R 

NAVA!,  2.  P 

ACE  SURVEILLANCE 

S  YSTEM 

PAHLCREN 

,  VA 

2  2  9  14  8 

01  0  Y 

ATTN 

C  A  P  T 

J.H.  BURTON 
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Or FICER-IN-CHARGE 
NAVAL  SURFACE  WEAPONS  CENTER 
WHITE  OAK,  SILVER  SPRING,  MD  2091 C 
0  1  C  Y  ATTN  CODE  FBI 

D I  RECTOR 

STRATEGIC  SYSTEMS  PROJECT  OFFICE 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  DC  20376 
0 1 C  Y  ATTN  NSP-21  A  1 

0 1 C  Y  ATTN  NSSP-2722  FRED  WIMBERLY 
COMMANDER 

NAVAL  SURFACE  WEAPONS  CENTER 
DAHLGREN  LABORATORY 
DAHLGREN,  VA  2  2  *4  *4  8 

0  1  C  Y  ATTN  CODE  D  F  -  1  R.  BUTLER 

OFFICER  OF  NAVAL  RESEARCH 
ARLINGTON,  VA  22217 
0 1 C  Y  ATTN  CODE  *465 
0 1 C  Y  ATTN  CODE  *4  61 
0  1  C  Y  ATTN  CODE  *402 
0 1  C  Y  ATTN  CODE  *4  2  0 
0  1  C  Y  ATTN  CODE  *421 


COMMANDER 

AEROSPACE  DEFENSE  COMMA  ND/ XPD 
DEPARTMENT  OF  THE  AIR  FORCE 
ENT  A  FB ,  CO  8091  2 
0  1  C  Y  ATTN  XPDQQ 


01  CY 

ATTN 

XP 

AIR  FORCE  GE 

OPHYSICS  LABORATORY 

HANSCOM 

A  FB  , 

MA 

01  731 

01  CY 

ATTN 

OPR 

HAROLD 

GARDNER 

0  1  C  Y 

ATTN 

L  KB 

KENNETH 

S.W.  CHAMPION 

0  1  C  Y 

ATTN 

OPR 

ALVA  T. 

STAIR 

01  CY 

ATTN 

PHD 

JURGEN 

BUCHAU 

0  1  C  Y 

ATTN 

PHD 

JOHN  P. 

MULLEN 

AF  WEAPONS  LABORATORY 

KIRTLAND 

AFT 

,  NM 

87117 

01  CY 

ATTN 

SUL 

01  CY 

ATTN 

CA 

ARTHUR  H 

.  GUENT1 

A  FTAC 

PATRICK 

A  FB  , 

FL 

32925 

0  1  C  Y 

ATTN 

TN 

AIR  FORCE  AVIONICS  LABORATORY 
WRIGHT  - PATTERSON  AFE,  OH  -8*433 
0 1 C  Y  ATT  N  A  AD  W ACE  HUNT 
C1CY  ATTN  A  AT  ALLEN  JOHN SC 

DEPUTY  CHIEF  OF  STAFF 
RESEARCH,  DEVELOPMENT,  5,  ACC 
DEPARTMENT  OF  THE  AIR  FORCE 
WASHINGTON,  DC  20330 
0  1  C  Y  ATTN  AFRDC 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  A  F  3 ,  MA  01731-  5000 
0 1 C  Y  ATTN  J.  DEAS 
ES  D  /  S  C  D  -  -4 

COMMAN  DER 

FOREIGN  TECHNOLOGY  DIVISION,  A 
WRI  GHT  - PATTERSON  A  F3  ,  OH  *4  5  *4  3  3 
0 1  C  Y  ATTN  N I  CD  LI3RARY 
0  1  C  Y  ATTN  ET  DP  3.  BALLARD 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  A 
GRIFFISS  A  F3 ,  NY  13*4*41 

01  C  Y  ATTN  DOC  LI3RARY/TSLD 
0  1  C  Y  ATTN  OCSE  V.  COYNE 

STRATEGIC  AIR  COMMAND/XPFS 
OFFUTT  A  FB ,  NB  681  1 3 
0 1 C  Y  ATTN  XPFS 

SAMSO/MN 

NORTON  A  F3  ,  CA  9  2  -*  C  9 
(MINUTEMAN) 

0 1  C  Y  ATTN  MNNL 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  A 
HANSCOM  A  F  3 ,  MA  01731 

0 1 C  Y  ATTN  EEP  A.  LORENTZEN 

DEPARTMENT  OF  ENERGY 
LIBRARY  ROOM  C-0*42 
WASHINGTON,  DC  205^5 

01  CY  ATTN  DOC  CON  FOR  A.  LA 

DEPARTMENT  OF  ENERGY 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.O.  BOX  5*400 
ALBUQUERQUE ,  NM  97 1 1 5 

01CY  ATTN  DOC  CON  FOR  D.  SH 


I 


EG&G ,  INC. 

LOS  ALAMOS  DIVISION 

P . 0 .  BOX  309 

LOS  ALAMOS.  NM  355  A U 

J  1C  Y  ATTN  DOC  CON  FOR  J.  BREEDLOVE 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  9^550 

0 1 C  Y  ATTN  DOC  CON  FOR  TECH  INFO 
DEPT 

0 1  C  Y  ATTN  DOC  CON  FOR  L-3S9  R.  OTT 
0 1 C  Y  ATTN  DOC  CON  FOR  L -  3 1  R.  HAGER 

LOS  ALAMOS  NATIONAL  LABORATORY 

P.O.  BOX  1663 

LCS  ALAMOS,  NM  875^5 


0  1  C  Y 

ATTN 

DOC 

CON 

FOR 

J  . 

WOLCOTT 

0  1  C  Y 

ATTN 

DOC 

CON 

FOR 

R  . 

F.  TASCHEK 

0  1  C  Y 

ATTN 

DOC 

CON 

FOR 

E. 

JONES 

01  CY 

ATTN 

DOC 

CON 

FOR 

J  . 

MALIK 

0  1  C  Y 

ATTN 

DOC 

CON 

FOR 

R. 

JEFFRIES 

01  C  Y 

ATTN 

DOC 

CON 

FOR 

J  . 

ZINN 

0  1  C  Y 

ATTN 

DOC 

CON 

FOR 

D  . 

WESTERVELT 

01  C  Y 

ATTN 

D. 

SAPPENFIELD 

LOS  ALAMOS  NATIONAL  LABORATORY 
MS  DH38 

LOS  ALAMOS,  NM  875^5 
0 1 C  Y  ATTN  S.P.  GARY 
0 1 C  Y  ATTN  J.  BOROVSKY 

SANDIA  LABORATORIES 
P.O.  BOX  5800 
ALBUQUERQUE,  NM  87115 

01CY  ATTN  DOC  CON  FOR  W.  BROWN 
0 1 C  Y  ATTN  DOC  CON  FOR  A. 

THORNBROUGH 

0 1 C  Y  ATTN  DOC  CON  FOR  T.  WRIGHT 

0 1  C  Y  ATTN  DOC  CON  FOR  D.  DAHLGREN 

0 1 C  Y  ATTN  DOC  CON  FOR  3 ! 41 

0 1 C  Y  ATTN  DOC  CON  FOR  SPACE  PROJECT 
DIV 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  94550 

0 1 C  Y  ATTN  DOC  CON  FOR  B.  MURPHEY 

0 1 C  Y  ATTN  DOC  CON  FOR  T.  COOK 


OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  DC  2054  5 

0 1 C  Y  ATTN  DOC  CON  DR.  YO  SONG 

NATIONAL  OCEANIC  &  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 

0 1 C  Y  ATTN  R.  GRU3B 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 


AEROSPACE  CORPORATION 


P.O.  B 

OX  92957 

LOS  AN 

GELES  , 

C  A 

90009 

01  CY 

ATTN 

I  . 

GARFUNKEL 

0  1  C  Y 

ATTN 

T  . 

SALMI 

01  CY 

ATTN 

V  . 

JOSEPHSON 

0  1  C  Y 

ATTN 

S. 

BOWER 

01  CY 

ATT  N 

D  . 

OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

0 1 C  Y  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOC.,  INC. 

1901  RUTLAND  DRIVE 
AUSTIN,  TX  78758 

01CY  ATTN  L.  SLOAN 
0 1 C  Y  ATTN  R.  THOMPSON 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  BOX  983 
BERKELEY,  CA  94701 

0 1 C  Y  ATTN  J.  WORKMAN 
0 1 C  Y  ATTN  C.  PRETTIE 
0 1 C  Y  ATTN  S.  BRECHT 

BOEING  COMPANY  ,  THE 
P.O.  BOX  3707 
SEATTLE,  WA  98124 


01C  Y 

ATTN 

G  . 

KEISTER 

01  CY 

ATTN 

D  . 

MURRAY 

0  1  C  Y 

ATTN 

G  . 

HALL 

01  CY 

ATTN 

J  . 

KENNEY 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE ,  MA  021 3? 

0 1 C  Y  ATTN  D.  B.  COX 
01C Y  ATTN  J . P .  GILMORE 


COMSAT  LABORATORIES 
22300  COMSAT  DRIVE 
CLARKSBURG,  MD  20371 
0  1  C  Y  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

0 1 C  Y  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1 359 

RICHARDSON,  TX  75080 
0 1 C  Y  ATTN  H.  LOGSTON 
0 1 C  Y  ATTN  SECURITY  (PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

606  Wilshire  Blvd. 

Santa  Monica,  CA  90401 
0 1 C  Y  ATTN  C.B.  GA3BARD 
0 1 C  Y  ATTN  R.  LELEVIER 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 
0 1 C  Y  ATTN  M.H.  BORTNER 
SPACE  SCI  LAB 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA' 

FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 

0 1 C  Y  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
0 1 C  Y  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

GTE  SYLVAN  IA,  INC  . 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM  ,  MA  021  9*1 

0 1 C  Y  ATTN  DICK  STEINHOF 

HSS ,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

0 1 C  Y  ATTN  DONALD  HANSEN 


ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61 820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 
01CY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA ,  VA  2231 1 
0 1 C  Y  ATTN  J.M.  A  E I N 
0 1  C  Y  ATTN  ERNEST  BAUER 
0  1  C  Y  ATTN  HANS  WOLFARD 
0  1  C  Y  ATTN  JOEL  BENGSTON 

I  NT  L  TEL  &  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY ,  NJ  07110 

0 1 C  Y  ATTN  TECHNICAL  LIBRARY 

JA  YCOR 

11011  TORRE  YANA  ROAD 
P.O.  BOX  8515^ 

SAN  DIEGO ,  CA  92138 

0 1 C  Y  ATTN  J.L.  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

0 1 C  Y  ATTN  DOCUMENT  LIBRARIAN 
0 1 C  Y  ATTN  THOMAS  POTEMRA 
01CY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
0 1 C  Y  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED 
STUDI ES 

816  STATE  STREET  (P.O  DRAWER  QQ  ) 
SANTA  BARBARA,  CA  93102 
0 1 C  Y  ATTN  DAS  I  A  C 
0 1 C  Y  ATTN  WARREN  S.  KNAPP 
0 1 C  Y  ATTN  WILLIAM  MCNAMARA 
0 1 C  Y  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 

0 1 C  Y  ATTN  IRWIN  JACOBS 


61 


LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
P.O.  BOX  50*4 
SUNNYVALE,  CA  9*4086 

0 1 C  Y  ATTN  DEPT  60-12 
0 1 C  Y  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC. 
3251  HANOVER  STREET 
PALO  ALTO,  CA  9*4  3  0*4 

0 1 C  Y  ATTN  MARTIN  WALT  DEPT  52-12 
0  1 C  Y  ATTN  W.L.  IMHOF  DEPT  52-12 
01CY  ATTN  RICHARD  G.  JOHNSON 
DEPT  52-12 

01CY  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 

0 1 C  Y  ATTN  R.  HEFFNER 

MCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  926*4  7 
0 1 C  Y  ATTN  N.  HARRIS 
01CY  ATTN  J.  MOULE 
0 1 C  Y  ATTN  GEORGE  MROZ 
0 1 C  Y  ATTN  W,  OLSON 
0 1 C  Y  ATTN  R.W.  HALPRIN 
0 1 C  Y  ATTN  TECHNICAL 

LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
0 1 C  Y  ATTN  P.  FISCHER 
0 1 C  Y  ATTN  W.F.  CREVIER 
0 1 C  Y  ATTN  STEVEN  L.  GUTSCHE 
0 1 C  Y  ATTN  R.  BOGUSCH 
0 1 C  Y  ATTN  R.  HENDRICK 
0 1 C  Y  ATTN  RALPH  KILB 
0 1 C  Y  ATTN  DAVE  SOWLE 
0 1 C  Y  ATTN  F.  FAJEN 
0 1 C  Y  ATTN  M.  SCHEIBE 
0 1 C  Y  ATTN  CONRAD  L.  LONGMIRE 
0 1 C  Y  ATTN  B.  WHITE 
0 1 C  Y  ATTN  R.  STAGAT 

MISSION  RESEARCH  CORP. 

1720  RANDOLPH  ROAD,  S.E. 

ALBUQUERQUE,  NM  87106 
0 1 C  Y  R.  STELLINGWERF 
0 1 C  Y  M.  ALME 
01CY  L.  WRIGHT 


MITF.E  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  V A  22101 

0 1 C  Y  ATTN  W.  HALL 
0 1 C  Y  ATTN  W.  FOSTER 

PACIFIC -SIERRA  RESEARCH  CORP 
1  23*40  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

01 C  Y  ATTN  E.C .  FIELD,  JR. 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

0 1 C  Y  ATTN  IONOSPHERIC  RESEARCH  LA 

PHOTOMETRICS,  INC. 

*4  ARROW  DRIVE 
WOBURN,  MA  C1801 

01 C  Y  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

0 1 C  Y  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  9*4610 
ATTN  A.  THOMSON 

R  &  D  ASSOCIATES 

P.O.  BOX  9695 

MARINA  DEL  REY,  CA  90291 

0 1 C  Y  ATTN  FORREST  GILMORE 
0 1 C  Y  ATTN  WILLIAM  B.  WRIGHT,  JR. 

0 1 C  Y  ATTN  WILLIAM  J.  KARZAS 
01CY  ATTN  H.  OR Y 
0 1 C  Y  ATTN  C.  MACDONALD 
0 1 C  Y  ATTN  BRIAN  LAMB 
01CY  ATTN  MORGAN  GROVER 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

01  CY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  *4  2  n  d  STREET 
NEW  YORK,  NY  10036 

01  GY  ATTN  VINCE  TRAPANI 


SCIENCE  APPLICATIONS 

INTERNATIONAL  INCORPORATED 


1  1  50  PROSPECT  PLAZA 


LA  JOLLA 

,  C  A 

92037 

01  C  Y 

ATTN 

LEWIS  M.  LI NS  ON 

0  1  C  Y 

ATTN 

DANIEL  A.  HAMLIN 

01  CY 

ATTN 

E.  FRIEMAN 

0  1  C  Y 

ATTN 

E  .  A .  STRAKER 

01  CY 

ATTN 

CURTIS  A.  SMITH 

SCIENCE  APPLICATIONS 

INTERNATIONAL  CORPORATION 
1 71 0  GOODRIDGE  DR. 

MCLEAN ,  VA  221 02 
0  1  C  Y  J.  COCKAYNE 
0 1 C  Y  E.  HYMAN 


SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 


MENLO 

P  ARK  , 

CA  94025 

01  CY 

ATTN 

J.  CASPER 

0  1  C  Y 

ATTN 

DONALD  NEILSON 

01  C  Y 

ATTN 

ALAN  BURNS 

0  1  C  Y 

ATTN 

G.  SMITH 

01  C  Y 

ATTN 

R.  TSUNODA 

0  1  C  Y 

ATTN 

DAVID  A.  JOHNSON 

01  CY 

ATTN 

WALTER  G.  CHESNUT 

0  1  C  Y 

ATTN 

CHARLES  L.  RINO 

01  CY 

ATTN 

WALTER  JAYE 

0  1  C  Y 

ATTN 

J.  VICKREY 

01  CY 

ATTN 

RAY  L.  LEADABRAND 

0  1  C  Y 

ATTN 

G.  CARPENTER 

01  CY 

ATTN 

G.  PRICE 

0  1  C  Y 

ATTN 

R.  LIVINGSTON 

01  CY 

ATTN 

V.  GONZALES 

0  1  C  Y 

ATTN 

D.  MCDANIEL 

V  I S  I  D  Y  N  E 

SOUTH  BEDFORD  STREET 
BURL INGTON ,  KA  01803 
0 1 C  Y  ATTN  W  .  REIDY 
0 1 C  Y  ATTN  J.  CARPENTER 
0 1  C  Y  ATTN  C.  HUMPHREY 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  1 521 3 

0 1 C  Y  ATTN:  N.  ZABUSKY 


Code  1220  lcp 
Records  lcp 

Director  of  Research 
U.S.  Naval  Academy 
Annapolis,  MI)  21402  2cp 


TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

0 1 C  Y  ATTN  W. P.  BO  QU 1ST 

TRW  DEFENSE  &  SPACE  SYS  GROUP 


ONE  SPACE  PARK 


REDONDO 

BEACH 

• 

CA  90278 

01CY 

ATTN 

R  . 

K.  PLEBUCH 

01  C  Y 

ATTN 

S. 

ALTSCHULER 

0  1  C  Y 

ATTN 

D  . 

DEE 

01  CY 

ATTN 

D/ 

STOCKWELL 

SNTF/1 575 

63 


V 


IONOSPHERIC  HODELING  DISTRIBUTION  LIST 
(UNCLASSIFIED  ONLY) 


PLEASE  DISTRIBUTE  ONE  COPY  TO  EACH 
OTHERWISE  NOTED) 


NAVAL  RESEARCH  LABORATORY 


WASHINGTON,  DC  20375 

DR. 

H. 

GURSKY  -  CODE  4100 

DR. 

J. 

M.  GOODMAN  -  CODE  4180 

DR. 

P. 

RODRIQUEZ  -  CODE  4750 

DR. 

P. 

MANGE  -  CODE  4101 

DR. 

R. 

MEIER  -  CODE  4140 

A. F.  GEOPHYSICS  LABORATORY 
L.G.  HANSCOM  FIELD 
BEDFORD,  MA  01731 
DR.  T.  ELKINS 
DR.  U.  SWIDER 
MRS.  R.  SAGALYN 
DR.  J.M.  FORBES 
DR.  T.J.  KENESHEA 
DR.  W.  BURKE 
DR.  H.  CARLSON 
DR.  J.  JASPERSE 
DR.  F.J.  RICH 
DR.  N.  MAYNARD 
DR.  D.N.  ANDERSON 
DR.  S.  BASU 

BOSTON  UNIVERSITY 
DEPARTMENT  OF  ASTRONOMY 
BOSTON,  MA  02215 
DR.  J.  AARONS 
DR.  M.  MENDILLO 

CORNELL  UNIVERSITY 
ITHACA,  NY  14850 
DR.  B.  FEJER 
DR.  R.  SUDAN 
DR.  D.  FARLEY 
DR.  M.  KELLEY 

INSTITUTE  FOR  DEFENSE  ANALYSIS 
1801  N.  BEAUREGARD  STREET 
ARLINGTON,  V A  22311 
DR.  E.  BAUER 


OF  THE  FOLLOWING  PEOPLE  (UNLESS 


MASSACHUSETTS  INSTITUTE  OF  TECHNO! OGY 
PLASMA  FUSION  CENTER 
CAMBRIDGE,  MA  02139 
LIBRARY,  NV16-262 
DR.  T.  CHANG 
DR.  R.  LINDZEN 

NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREENBELT,  MD  20771 

DR.  N.  MAYNARD  (CODE  696) 

DR.  R.F.  BENSON 
DR.  K.  MAEDA 
DR.  S.  CURTIS 
DR.  M.  DUBIN 

COMMANDER 

Naval  ocean  systems  center 
SAN  DIEGO,  CA  92152 

MR.  R.  ROSE  -  CODE  5321 

NOAA 

DIRECTOR  OF  SPACE  AND 
ENVIRONMENTAL  LABORATORY 
BOULDER,  CO  80302 
DR.  A.  GLENN  JEAN 
DR.  G.V.  ADAMS 
DR.  K.  DAVIES 
DR.  R.F.  DONNELLY 

OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VA  22217 
DR.  G.  JOINER 

LABORATORY  FOR  PLASMA  AND 
FUSION  ENERGIES  STUDIES 
UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20742 
JHAN  VARYAN  HELLMAN, 

REFERENCE  LIBRARIAN 


PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PA  16802 


DR. 

DR. 


J.S.  NISBET 
P.R.  ROHRBAUGH 


PRINCETON  UNIVERSITY 
PLASMA  PHYSICS  LABORATORY 
PRINCETON,  NJ  08540 
DR.  F.  PERKINS 


SAIC 

1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 
DR.  D. A.  HAMLIN 
DR.  L.  LINSON 
DR.  E.  FRIEMAN 


SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 
DR.  R.  TSUNODA 

WALTER  CHESNUT 


DR. 


UNIVERSITY  OF  CALIFORNIA 
LOS  ALAMOS  NATIONAL  LABORATORY 
EES  DIVISION 
LOS  ALAMOS,  NM  87545 


DR. 

L. 

A.  CARPENTER 

DR. 

M. 

PONGRATZ,  ESS-DOT 

DR. 

M. 

LEE 

DR.  D. 

SIMONS,  ESS-7,  MS 

DR. 

R. 

DIVANY 

DR. 

L. 

DUNCAN,  ESS-7,  MS 

DR. 

P. 

BENNETT 

DR.  P. 

BERNHARDT,  ESS-7, 

DR. 

E. 

BLEVANS 

DR. 

S. 

P.  GARY,  ESS-8 

DENNIS 

RIGGIN,  ATMOS  SCI 

UNIVERSITY  OF  ILLINOIS 
DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
1406  W.  GREEN  STREET 
URBANA,  IL  61801 
DR.  ERHAN  KUDEKI 


UNIVERSITY  OF  CALIFORNIA, 
LOS  ANGELES 
405  HILLGARD  AVENUE 
LOS  ANGELES,  CA  90024 
DR.  F.V.  CORONITI 
DR.  C.  KENNEL 
DR.  A. Y.  WONG 


UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20740 


DR.  CHARLES  RINO 

DR.  K.  PAPADOPOULOS 

DR.  J.  VICKREY 

DR.  R.  LIVINGSTON 

DR.  E.  OTT 

JOHNS  HOPKINS  UNIVERSITY 

STANFORD  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

STANFORD,  CA  94305 

JOHNS  HOPKINS  ROAD 

DR.  P.M.  BANKS 

LAUREL,  MD  20810 

DR.  R.  HELLIUELL 

DR.  R.  GREENVALD 

DR.  C.  MENG 

U.S.  ARMY  ABERDEEN  RESEARCH 

AND  DEVELOPMENT  CENTER 

DR.  T.  POTEMRA 

BALLISTIC  RESEARCH  LABORATORY 

UNIVERSITY  OF  PITTSBURGH 

ABERDEEN,  MD 

PITTSBURGH,  PA  15213 

DR.  J.  HEIMERL 

DR.  N.  ZABUSKY 

DR.  M.  BIONDI 

GEOPHYSICAL  INSTITUTE 

UNIVERSITY  OF  ALASKA 

DR.  E.  OVERMAN 

FAIRBANKS,  AL  99701 

UNIVERSITY  OF  TEXAS  AT  DALLAS 

DR.  L.C.  LEE 

CENTER  FOR  SPACE  SCIENCES 

P.O.  BOX  688 

UTAH  STATE  UNIVERSITY 

RICHARDSON,  TX  75080 

4TH  AND  8TH  STREETS 

DR.  R.  HEELIS 

LOGAN,  UT  84322 

DR.  V.  HANSON 

DR.  R.  HARRIS 

DR.  K.  BAKER 

DR.  R.  SCHUNK 

DR.  J.  ST. -MAURICE 

DR.  N.  SINGH 

DR.  J.P.  McCLURE 
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